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1. Introduction 



Lanchester-type* differential-equation combat models are an important 
tool for analyzing many important problems of military operations research. 

In such a combat model, a so-called attrition-rate coefficient represents 
the fire effectiveness of a particular weapon-system type against a particular 
target type, i.e. the weapon-system type’s effective firepower against such 
a target. Time-dependent attrition-rate coefficients are used to model 
temporal variations in firepower on the battlefield. Thus, we see that 
time-dependent attrition-rate coefficients are important (and, in fact, 
essential [4-6]) for the quantitative analysis of hypoLhet i cal combat. 

Militarily realistic computer-based Lanchester-type models of quite 
complex military systems have been developed for almost the entire spectrum 
of combat operations, from combat between battalion-sized units to theater- 
level operations. Nevertheless, a simple combat model may yield a clearer 
understanding of significant interrelationships that are difficult to per- 
ceive in a more complex model, and such insights can subsequently provide 
valuable guidance for more detailed computerized investigations. In this 
report we consider such a simplified variable-coefficient Lanchester-type 
model of combat between two homogeneous forces. 

For this variable-coefficient Lanchester-type model of combat 
between two homogeneous forces, different functional forms for the attrition- 
rate coefficients lead to different mathematical functions being involved 
in representing and computing the force-level trajectories. In a previous 
paper [5] we have discussed the plausibility of the hypothesis that except 
for the special case of a constant ratio of attrition-rate coefficients, 

* 

So-called after pioneering work of F. W. Lanchester [3]. 
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the solutions to such differential equations cannot be represented in term 
of "elementary" functions of analysis. Thus, new transcendental functions 
arise in the study of combat modelled with time-dependent attrition-rate 
coefficients. In particular, we have previously introduced [5-6] so-called 
Lanchester-Clif ford-Schlaf li (LCS) functions for analyzing combat modelled 
with power attrition-rate coefficients with "no offset" (see Section 3 below). 

In the Appendix to this report is contained a reduced set of tables 
for the LCS functions: it contains tables of five-decimal-place values of 

the hyperbolic-like LCS functions F (x) , H. (x) , and T (x) (see Section A 

a 1-a a 

below) for 11 fractional values of a (see Section 6 below) . A companion 
report [8] contains the most extensive set of tables currently available. 

The main body of this report provides the theoretical and modelling background 
for the use of these tables. In particular, we examine a model of a constant- 
speed attack on a static defensive position and show how associated range- 
dependent kill rates give rise to time-dependent attrition-rate coefficients 
with "no offset." Numerical computations are presented to illustrate the 
use of the LCS functions for analyzing such "aimed-fire" combat. As a con- 
sequence of the availability of these tables, one can now study this variable- 
coefficient combat model almost as easily and thoroughly as Lanchester ! s 
classic constant-coefficient model. 



2 . Variable-Coefficient Lanches ter-Type Equations of Modern Warfare . 

We consider combat between two homogeneous forces modelled by the 
following var iable-coef f icient Lanches ter- type [3] (see [A, 5]) equations 
of modern warfare 
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with x(0) = Xq , 



( 2 . 1 ) 




with y (0) = y Q , 



where t = 0 denotes the time at which the battle begins, x(t) and y(t) 
denote the numbers of X and Y at time t, and a(t) and b(t) denote 
time-dependent Lanchester attrition-rate coefficients, which represent the 
effectiveness of each side*s fire. These coefficients depend on variables 
such as force separation, tactical posture of targets, rate of target acqui- 
sition, firing rate, etc. (see [4-7] for further details). Variable attrition- 
rate coefficients are used to model temporal variations in firepower on the 
battlefield. In any analysis of combat, moreover, we should use the above 
equations (2.1) only for x and y > 0 and, for example, set dx/dt = 0 
when x = 0, since negative force levels have no physical meaning. 



Mathematically, we assume that the attrition-rate coefficients a(t) 



and b(t) are defined, positive, and continuous for t^ < t < + °° with 
t 0 — We a ^ so assume that a(t) and b(t) € L(t^,T) for any finite 

T >_ tp. We further take a(t) and b(t) to be given in the form 



a (t) = k g(t) , 
a 



and 



b (t ) = k b h(t) , 



( 2 . 2 ) 



where k 




a 



when g(t) r h(t) . We introduce the combat-intensity parameter A^ and 

the relative-fire-effectiveness parameter A defined by 

R 
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From our assumptions about a(t) and b(t), it follows that, for example, 

T 

a(t) f L(t ,T) implies / a(t)dt = + <» . 

u t o 

The X force level as a function of time may be represented 

as [5,6] 

x(t) = x Q {C Y (0)C x (t) - S y (0)S x (t)} - y Q /^ (C x (0)S x (t) - S x (0)C x (t)}, (2.4) 



where the hypeibolic-like general Lanchester functions (GLF) C (t) and S (t) 

X A 

are linearly-independent solutions to the X force-level equation 




1 da 
a(t) dt 



dx 

dt 



a (t)b (t)x 



0 , 



with initial conditions 



(2.5) 



c x (t o ) 1 ’ W ° ’ 

( 2 . 6 ) 

(l/a(t Q ) } dC x /dt(t Q ) = 0 , { 1/ a (t Q ) } dS x /dt(t 0 ) = 1//^ . 

Here t^ denotes the largest finite time at which a(t) or b(t) ceases 
to be defined, positive, or continuous. The Y force level as a function 
of time is given by a similar expression, with C^(t) and S^(t) being 
analogously defined for the corresponding Y force-level equation. 

It is sometimes convenient to introduce the new independent variable 



t defined by 



t = / /a (s)b (s) ds . 

t 0 



(2.7) 



It is readily seen that the transformation t = x(t) is well defined and 
invertible. Let us denote x(0) as x^. We observe that tg <_ 0 implies 
that tg _> 0. If we denote the "average intensity of combat" as 
/a (t)b (t) , then 



/a(t)b(t) t = { ( 1 / 1 ) | l/ a (s)b (s) ds}t = t- t 

0 



0 



( 2 . 8 ) 



The substitution (2.7) transforms (2.5) into 



A 

dx 2 



- <!>{£ 



•• x = 0 , 



(2.9) 



with initial conditions 



x(t q ) = x Q , and { 1/ /r( t q ) } dx/di^) = -y , 



where R(x) = a(t)/b(t) 



3. Combat Modelled with Power Attrition-Rate Coefficients . 

The above equations (2.1) basically apply to M aimed-fire n combat 
when target-acquisition times do not depend on the numbers of targets available 
(see [5,6] for further details). A large class of tactical situations of 
interest can be modelled with the following general power attrition-rate 
coefficients [5-7] 
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and 



b(t) = k^ (t + C + A) V , 



(3.1) 



a(t) = k (t + C) P 
a 



where A and C _> 0. We will call A the offset parameter , since it allows 
us to model (with p and v _> 0) battles between opposing weapon systems 
with different maximum effective ranges (see [5,6]). We will call C the 
starting parameter , since it allows us to model (again, with p and v >_ 0) 
battles that begin within the maximum effective ranges of the two opposing 
systems. We observe that for the general power attrition-rate coefficients 
(3.1) we have t^ = -C, and p and v must be > -1 in order that a(t) 
and b (t) £ L(t^ ,T) . 

The above nomenclature is motivated and possible applications of our 
work are indicated by considering S. Bonder’s model of the constant-speed 
attack on a static defensive position (see [4-7] for further details) 



77 = -a(r)y, and 7 ^ = -B(r)x , (3.2) 

at at 

where r denotes the range between opposing forces, and a(r) and 0(r) 
denote range-dependent attrition-rate coefficients. Range is related to time 
by 

r (t) = R Q - vt , (3.3) 

where R^ denotes the opening range of battle and v > 0 denotes the 
constant attack speed. For example, let us consider the constant-speed attack 
of a homogeneous Y force against the static defensive position of a homo- 
geneous X force. Figure 1 diagrammatically portrays this situation. 
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\ 
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Figure 1. Diagram of Bonder's constant-speed attack model. 

Force separation, r(t), is given by r(t) = R -vt. 



The basic idea is that force separation, i.e. range between the opposing 
forces, changes over time, and the fire effectiveness of, for example, 
a single Y firer, denoted as a(r), depends on this force separation. 

In many cases of tactical interest, we may model the fire effective- 
ness of, for example, the Y weapon system (as a function of range) with 



a(r) 




0 



for 0 < r < R , 
— — a 



for R < r , 
a — 



(3.4) 



where R^ denotes the maximum effective range of the Y weapon system and 

y _> 0. Here p is used to model the range dependency of Y f s attrition-rate 
coefficient (see Figure 2). We model 3(r) similarly, with corresponding 
quantities R and v being analogous to R^ and p above. 

If we use (3.3) to eliminate range r from (3.4), we obtain 



dx 

dt 



= -a (t) y , 



= ~ Mt)x ’ 



(3.5) 



where the time-dependent attrition-rate coefficients a(t) and b(t) are 
given by (3.1). It follows that the offset and starting parameters are given 
•by 
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A = 




and 



C 




(3.6) 



and that 




and = B q 




V 



(3.7) 




ing (3.6) and Figure 3, the reader should have no trouble in understanding 
our terminology for A and C. 

When the offset parameter is equal to zero (i.c. A = 0) , then the 
coefficients (3.1) reduce to 



We will refer to (3.8) as power attrition-rate coefficients with "no offset . 11 
As we have seen above in Bonder’s constant-speed attack model, these 
coefficients model, for example, combat between weapon systems with the same 
maximum effective range so that there is no "offset 11 in the "reaching out" 
of the weapon systems against each other in combat (again, see Figure 3). 

For these coefficients (3.8), the transformed X force-level equation (2.9) 
becomes 



with initial conditions 



a(t) = k (t + C) y , 
a 



and 



b ( t) = (t + C) V . 



(3.8) 
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X ( T 0> = x o ’ and 



l(I ) 2q_1 

V dT 

[ T = T . 



2q-l 



= -y„/xT" 



0 R \ p + v + 2 



Here 



q = 



v + 1 
U + v + 2 



(3.10) 



T 



' 2A I 

M + v + 2 



(t + C) 



(U+ v+2)/2 



(3.11) 



and 



_ I 2X I \ (m+v+2)/2 

T 0 ^ p + v + 2 J 



(3.12) 



Let us observe that 0 < q < 1 when p and v > -1 . Furthermore, q > 1/2 
if and only if dR/dt < 0, i.e. R(t) is a strictly decreasing function of 
time . 



4 . Lanchester-Clif f ord-Schlaf li (LCS) Functions . 

Consider the function F^ (x) defined by the power series 



V x) 



r(a) 



y (x/ 2) 2k 

{k! r (k + a) } 



(4.1) 



For a ^ 0, -1, -2,... the radius of convergence for F a (x) is infinite 
by the ratio test for convergence of power series [2], Hence, F a (z) is 
an entire function of the complex variable z = x + iy, with an essential 
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singularity at the point at infinity. Now consider the function H^(x) 
defined by the infinite series 



H (x) 

a 



r(a) j 

k=0 



(x /2) 2(k+a) 

{k! r (k + a +1)} ' 



(4.2) 



Observing that 



H (x) = (1/a) (x/2) 2a F (x) , 
a a+1 



(4.3) 



we see that for a > 0 the infinite series (4.2) is uniformly convergent 
on compact subsets of the complex plane. From (4.3) one can readily deduce 
the recursive relation 



F 

a 



(x) 



F « + l< x) 



(x/2) 2 I 

| ot(a + 1) j 



f a + Z M 



(4.4) 



We will call the functions F (x) and H (x) Lanchester-Cli f ford-Schlaf li 

a a 

(LCS) functions (see Note 10 on pp . 66-67 of [5]). Other properties are 
readily deduced and are given in Table I, 

The function F^(x) satisfies the linear second-order ordinary 
differential equation 



d 2 F . . dF 

a .2a-l a 

O ' v / 



dx 



dx 



- F = 0 , 
a 



with initial conditions 



(4.5) 
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Table I 



1 . 



2 . 



3. 



4. 



5. 



6 . 



7. 



8 . 



9. 



Properties of the LCS Functions F (x) and H (x) . 

a a 



dF /dx = (x/2)^ ^ a H (x) 
a a 



dH /dx = (x/2)^ 3 ^F (x) 
a a 



F (x)F (x) -H (x)H 1 (x) = 1 Vx 

a v ' 1-a a' 1-a 

where a is not an integer (including zero) 



F q (x=0) = 1 



H (x=0) = 0 for a > 0 
a 



dF a /dx(x=0) = 0 



{(x/2)^ ^ a dH /dx) n = 1 
a x=0 



Fp/ 2 (x) = cosh x 



H 1/2 (x) 



s ir.h x 
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v°> ■ 



and 



dF 
c 

dx 



( 0 ) = 0 , 



while H (x) satisfies 

a 



d 2 H , , dH 

- (*£i) - „ . 0 , 

,2 x dx a 

dx 



(4.6) 



with initial conditions 



H a (0) = 0 , 



and (~) 



( 1-2 a dH ) 



dx 



) x=0 



= 1 



Thus, {F , H, } is a fundamental system of solutions to 
a 1-a 



M + (— ) £ - * - 0 , 

,2 x dx 

dx 



(A. 7) 



with Wronskian W(F , IL ) = (x/2)^ ^ a , It follows that the GLF for the 

a 1-a 

X and Y force-level equations for combat modelled with the attrition-rate 
coefficients (3.8) are given by 



c x (t) ■ V T<t) )* s x (t) ■( t + V + t ) 2 ' 1 " 1 V T<t)) • «- 8 > 



C Y^ F p(x(t))» S Y^ \u + v + 2 



\l-2q 

) V* 



(t)) , (A. 9) 



where p = 1-q. If we define 
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(4.10) 



V x) = H l-a (x)/ V x) 



then 



S (t) / X \2q-l H (t (t) ) 

T x (t) = cTTtT = \ V + V + 2 ) /(t(0) ’ 

A q 



(4.11) 



or 



T X U> 



\i + v + 2 



\2q-l 

j ’ 



(4.12) 



where T (t) denotes a hyperbolic-like GLF, which corresponds to the 
A 

hyperbolic tangent. Observing that for p, v > -1, lim + ^ x(t) = + 
we see that T^(x) is a strictly increasing function of x on the interval 
[0, + ») and 



0 < T (x) 
— a 



r(l-q) 

r(a) 



for 0 < x < + °° , 



(4.13) 



with 



lira T (x) 
a 

x -y + co 



r(l-a) 
r(a) ’ 



(4.14) 



since by the results of Taylor and Comstock [7] the parity-condition parameter 
Q* = Q*(p, v, C = 0) is given by 



lira 
t-*- + 00 



T x (t) 



1 

Q*(p, v, 0) 




(4.15) 



We recall that Taylor and Comstock [7] have introduced the so- 
called parity-condition parameter Q* as the value (or range of such values) 
for the initial condition Q to the initial-value problem 
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with E x (t 0 ) = ] , 



(4.16) 



with E Y (t Q ) = Q , 

such that E (t;Q*) and E (t;Q*) > 0 for all t _> t . In other words, 

A I U 

Q* is the value of Q in (4.16) above such that neither E nor E 

A i 

ever become zero. In this case, both E (t;Q*) and E (t;Q*) are positive, 
strictly decreasing functions, similar to decreasing exponentials. Thus, 
we may call Q* "the Y equivalent of an X force of unit strength," 
since the forces are "at parity," with neither force being annihilated in 
finite time. Taylor and Comstock have shown that for either a(t) f. L(0,+ 
or b(t) f. L(0, + «>) , then Q* is unique and given by 



dE 

3T * - ^ 



dE Y r~ 

JC= - ^ b(t)E x 



S X (t) i 

i + . Cx(t) ’ Q* 



(4.17) 



The significance of the parity-condition parameter Q* is that it allows 
us to predict force annihilation as the following theorem shows. 



THEOREM 1 (Taylor and Comstock [7]): Assume that either 

a(t) £ L(0, + 00 ) or b(t) t L(0, + 00 ) . Then the X force 
will be annihilated in finite time if and only if 





c x (0) - 0*S X (0) | 
Q*C y (0) - S y (0) j 



(4.18) 
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5. Use of LCS Functions for Analyzing Combat . 



The Lanchester-Clif ford-Schlaf li (LCS) functions F (x) and H (x) 

a a 

are useful for analyzing ”aimed-f ire” combat (see Section 3 above) modelled 
with the power attrition-rate coefficients with "no offset” (3.8), which 
we rewrite here as 



a(t) = k^(t + C) U , and b(t) = k^t + C) V . 



(5.1) 



In other words, the LCS functions arise in solving the differential combat 
model (2.1) with at trit ion -rate coefficients (5.1). In order that both 
a(t) and b(t) € L(t^,T), we must have p and v > -1. Military situations 
modelled by these equations have been discussed in Section 3 above, e.g. 
combat between two weapon systems with the same maximum effective range. 

For such combat, the LCS functions may be used to 

(1) compute force-level declines, 

(2) predict force annihilation, 

and (3) predict the time of force annihilation. 

Let us now see how the LCS functions may be used to obtain the 
above information about force-level declines and force-annihilation pre- 
diction. According to (2. A), (A. 8), and (A. 9) above, the X force level 
is given by 



x(t) = x {F (t ) F (t (t) ) - H (t ) II (t (t ) ) } 
0 p 0 q q 0 p 



- y 



0 ^ [ v + v X + 2 / q {F q (T 0 )H p (T(t))_1I p (T 0 )F q (T(t))1 ’ (5 ‘ 2) 
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where q is given by (3.10), p = l~q, and 
we rewrite as 



T (t) is given by (3.11), which 



T ( t ) 



^ A J 

M ’+ v + 2 



(t 



+ c) (p+v+2) / 2 



(5.3) 



The time to annihilate the X force* is determined by x(t ) = 0, and 

3 

it follows that 



Xq (T (t j, , - 

WV + y 0 A rL + v\- 2, 



q-p 



iq-p 



H (T f ) 

_E_<L 

F (t ) 

q o 



(5.4) 



where from (4.10) 



T (i(t)) = 11 (t ( t ) ) /F ( t ( t ) ) , 

q p q 



(5.5) 



and we recall that p -f q = 1 . It follows that the time to annihilate X, 
X J 

t , is given by 
a 



If we multiply the first equation of (2.1) by y, the second by x, add, 
and integrate, we obtain 

t 9 9 

x ( t ) y (t) = x Q y Q - J {a (s) y (s) + b(s) x (s)}ds , 



which shows that x(t) and y(t) can have at most one finite zero. Hence, 
if x(t^) = 0, then we know that y(t) > 0 for all t >_ 0. 

3 
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(5.6) 



X -1 1 _-l 
t = t < T 
a \ q 



-oW+yo-'V (— ■ v V T ) q " H p ( V 



q-p 



* 0 yy + y o ^ (tt yv 



Taylor and Comstock [7] have shown that T (t) is strictly increasing and 

q 

satisfies (see also (A. 12) above) 



0 < T (t) < r(p)/r(q) , 
— q 



(5.7) 



where p = 1-q. It follows that in order for X to be annihilated in finite 
time, the right-hand side of (5.4) must be less than T(p)/r(q). Let us 
observe that for t = -C ■ 0, (5.4) simplifies to 



p-q 



T (x(O) = 

q a 



y o •y 



p 4- v + 2 



(5.8) 



Thus, we have proved the following theorem concerning force- 
annihilation prediction. 



THEOREM 2: Consider combat between two homogeneous forces 

modelled by (2.1) with pov/er attrition-rate coefficients 
(5.1). Assume that p and v > -1 and that the above 
equations hold for all time. Then the X force will be 
annihilated in finite time if and only if 
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r ( V 

R \ u + v + 2 / 



j 



r <<!) , * 0 F p cy + y 0 /' R \ --rvr-2 / VV| 



I 



‘ r(p) | x 0 H q <T 0 ) + y 0 \;r+TT7/ W I • 



q-p 



(5.9) 



where q = ( v + 1) / (p + v + 2) and p = 1-q. For t^ = 0 
(i.e. C = 0 so that Tq = 0) * X will be annihilated in finite 
time if and only if 



o < m /r ( — ^ — ) 

Yq r(q) R \ p + v + 2 / 



q-p 



(5.10) 



6 . Tabulation of LCS Functions . 

This report contains a reduced set of tables of the Lanchester-Clif ford- 
Schlafli functions. The Appendix contains tables of five-decimal-place values 

of the hyperbolic-like LCS functions F (x) , H, (x) , and T (x) for various 

a 1-a a 

values of the argument x, namely x = 0.00 (0.01) 2.00 (0.1) 10.0, and 
a = 1/2, 1/3, 2/3, 1/4, 3/4, 1/5, 2/5, 3/5, 4/5, 3/7, and 4/7. These values 
of the index a correspond to u, v = 0, 1, 2, and 3 in (3.8) and allow one 
to analyze, for example, a basic spectrum of range capabilities for weapon 
systems in the constant-speed-attack model of Section 3. These tables have 
been calculated by the recursive means given in Section 8 of [5]. A more 
extensive tabulation (namely, for a = 1/2, 1/3, 2/3, 1/4, 3/4, 1/5, 2/5, 

3/5, 4/5, 2/7, 3/7, 4/7, 5/7, 4/9, 5/9, 3/11, 5/11, 6/11, 8/11, 5/13, 8/13, 

5/17, 12/17, 5/21, and 16/21 corresponding to p , = 0, 1/4, 1/2, 1, 1 j, 2, 3) 
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is to be found in a companion report [8]. This companion report contains 
the most extensive set of tables of the Lanchester-Clif f ord-Schlaf li 
functions currently available. 

A representative tabulation of the hyperbolic-like LCS functions 
F^(x), ^(x), anc * T^(x) as gi ven i n > f° r example, Tables 8A and 8B of 

the Appendix for a = 3/5. The values of the argument x are the same as those 
used for the tabulation of the hyperbolic functions by Abramowitz and Stegun 
[1]. We observe from Table 8B and (4.13) that the limiting value of 
T^(x) as x + 00 (here a = 3/5) is quickly reached, with three-decimal- 
place accuracy already attained for x = 4.5. Moreover, the use of these 
tables (specifically, Tables 8A and 8B of the Appendix) for combat analysis 
is illustrated in the next section. 



7 . Numerical Examples 

In this section we examine a couple of numerical examples to show 
some of the insights that may be gained into the dynamics of combat between 
two homogeneous forces from our results (see also [6]). These examples 
illustrate the use of the LCS functions F^(x), H^_^(x), and T^(x) f° r 
analyzing "aimed-f ire" combat modelled with the power attrition-rate coefficients 
with "no offset" (5.1). As in [4-7], we consider S. Bonder f s model (3.2) 
for the constant-speed attack against a static defensive position. We will 
focus on the use of the LCS functions for predicting force annihilation, 
since the computing of force-level trajectories with Lanchester functions 
is adequately handled elsewhere (see [4-5]), 
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Let us accordingly consider the constant-speed attack of a homogeneous 

Y force against the static defensive position of a homogeneous X force 
(see Section 3 above for further modelling details, especially Figure 1). 

For our numerical computations, we assume that the fire effectiveness of the 

Y weapon system varies linearly with range, i.e. 



cf(r) = 




for 0 < r < R , 
— — a 



fcr R < r , 
a — ’ 



(7.1) 



and that the fire effectiveness of the X weapon system varies quadratically 
with range, i.e. 



e(r) 




for 0 < r < R , 

p 

for R g <_ r , 



(7.2) 



with R^ = R^, i.e. both weapon systems have the same maximum effective range. 
In other words, p = 1 in (3.4) and v = 2 for B(r). We consider a battle 
modelled by the input data given in Table II. In terms of time as the 
independent variable, the attrition-rate coefficients (7.1) and (7.2) become 
via (3.3) 



a(t) = k (t + C) 
a 



and 



b(t) = k b (t + C) 2 , 



(7.3) 
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Table II. Input Data for Numerical Examples 



Vi = 1, v = 2 

c*q = 0.06 X casualties/minute/Y firer 



= 0.6 Y casualties/minute/X firer 



R = R_ = 2000 meters 
a 6 



v = 5 miles/hour 



where R 



a ^ B* 



C = 



R - R n 
a 0 



k » 
a 



v 



and 



k b ' B 0( r. 



(7.4) 



From the input data given in Table II, we compute the parameter values shown 

in Table III, since the transformed X force-level equation is given by (3.9) 

with q = (v + l)/(u + v + 2) , p = 1-q, p = 1 , and v = 2. Thus, the X force 

level may be computed with F (x) and H, (x) with a = q = 3/5. Force- 

a 1— a 

annihilation prediction involves the limiting value of T (x) = H (t)/F (t) 

•j. 1-a a 

as T ^+ co . From Table 8B of the Appendix and Table III, we note the predicted 

agreement between r(l-cx)/r(a) and the limiting value of T (x) as x ^ + 00 

a 

[recall (4.13)] for a = q = 3/5. We now consider two cases: (I) = 2000 

meters, and (II) R^ = 1250 meters. 

When Rq = 2000 meters (see Figure 3 of [4]), we have C = 0 and 

= 0. The maximum time that the battle can last is t = R^/v = 14.91 minutes, 
0 max 0 

since at this time the attackers reach their final objective, i.e. the 
defender’s position (again, see Figure 1). We now consider the qualitative 
behavior of the u = 1, v = 2 force-level trajectory shown in Figure 3 of 
[4]. Theorem 2 tells us that the X force can be annihilated in finite time 
if and only if 



*0 < He! 

y o r(q) 

where q = 3/5 and p = 1-q. 
compute from (7.3) that the X 
and only if 



^ ( p + V + 2 / * ( 

Using the numerical values in Table III, we 
force can be annihilated in finite time if 
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Table III. Parameter Values for Numerical Examples 



k = 4.0233 x 10 3 X casualties/ (minute)^/Y firer 
a 

k, = 2.6979 x 10 Y casualties/ (minute) V /X firer 
b 

p = 2/5, q = 3/5 

r (p)/r (q) = 1.48951 
A = 0 
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When the X force can be annihilated, its annihilation time is given by 



(5,8), which we rewrite here as 



where 



T (i(t )) = 
q a 



> 0 ^ 



}j 4 v 4 2 



p-q 



T (t) 



I 2A 1 \ (|i+v+2)/2 

\ U + v + 2 f 



(7.5) 



(7.6) 



Thus, for the numerical values given in Table III, the time of annihilation 
of the X force is given by 



T (i(t X )) 
q a 



3.544 




(7.7) 



with q = 3/5. We will now illustrate further 

= 30. From (7.4) we see that the X force 

X 

time (but we must verify that t < t ) . In 

a — max 



computations for = 10 and 

can be annihilated in finite 
this case (7.7) becomes 



T (x (t X ) ) = 1.18122 . (7.8) 

q a 



x 

We must now determine i(t ) such that 

a 

polation methods and Tables 8A and 8B. 



T(t X ) = T hi. 18122) 
a q 

From Table 8A, we find 



by using inter- 



T (t) = 1.18172 

q 


for 


= 1.01 


T (t) = 1.17630 


for 


= 1.00 



q 



so that using linear interpolation, we obtain 

i(t*) = 1.009 , (7.9) 

cl 



whence use of (7.6) yields 

t^ = 14.24 minutes, (7.10) 

which is less than t = 14.91 minutes so that the defending X force is 

max 

indeed annihilated before the attacking Y force reaches its final objective. 
Since r(t) = R^ - vt, we find that force separation at the instant of 
annihilation of the X force is 

r^ = 89.8 meters . (7.11) 

a 

Further results may be computed in a similar fashion and are given in 
Table IV. 

When Rq = 1250 meters (see Figure 3 of [5]), we have C = 5.5923 

minutes, t_ = 0.0975, and t = R_/v = 9.32 minutes. In this case 
0 max 0 

Theorem 2 tells us that the X force can be annihilated in finite time 
if and only if 
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Table IV. Annihilation of the X Force as a Function 



of the Initial Force Ratio for = 2000 meters 



w 


t^ (minutes) 
a 


r^ (meters) 
a 


0.333 


14.24 


89.8 


0.250 


11.61 


443.2 


0.200 


10.19 


633.2 
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(7.12) 



x 



0 



y 0 



< 




h \ q - p im 

V + v + 2 / r (q) 




) - 



IM 

t(p) 



) - 



'(p) 

r(q) 
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with q = 3/5 and p = 1-q. Using linear interpolation, we obtain from 
Tables 7A and 8A of the Appendix that for the numerical values of Table III 



F (r ) = 1.006 , H (t ) = 0.044 , 

p u q u 

(t q ) = 1.004 , H p ( V = 0,223 » 



(7.13) 



so that (7.12) says that the X force can be annihilated if and only if 



0 

— < 0.382 . 

y o 



(7.14) 



When the X force can be annihilated, its annihilation time is given by 
(5.4), which we rewrite here as 



0 



T (t(O) = 
q a 



y o \ 



/ A \ p_q 

( — 1 — ) 
y p + v + 2 / 



F (O + H (t ) 

p 0 p 0 



p-q 



w + 






, y + v + 2 



W 



(7.15) 



whence for the data of Table III 



3.565u + 0.223 

T a (T(t: a )) = 0.156u_ + 1.004 ’ 



(7.16) 
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where = x^/y^. Let us also record here that (3.11) yields 



( {h + V + 2 } T \ 
v 2A I ) 



2/ (}i+v+2) 



- C . 



(7.17) 



We will again illustrate further computations for = 10 and = 30. 

From (7.14) we see that the X force can be annihilated in finite time (but 

X 

again we must investigate whether or not t <_ t ) . In this case (7.16) 

a max 

becomes 



T (t ( t X ) ) = 1.33651 , 
q a 



(7.18) 



whence Table 8A of the Appendix and linear interpolation yield 



T (0 = 1.397 , 
a 



(7.19) 



so that by (7.17) 



= 10.63 minutes . (7.20) 

a 

y 

Since t = R^/v = 9.32 minutes < t , we see that the attacking Y force 
max 0 a ’ 

overruns the defender ! s position before annihilation of the X force occurs. 

Thus, the battle ends with = x(t r ) > 0 and y r > 0 at t r = t = 

f f J f f max 

9.32 minutes. Corresponding to t^ = 9.32 minutes is = 1.1318, and 
then Table 8A of the Appendix yields 
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F (t. = 1.1318) = 1.589 , 
9 f 



H (1.1318) = 1.973 



(7.21) 



P 



whence via (2.4), (4.8), (4.9), and (7.13) we obtain 



x^. = x(tp = x(r = 0) = 1.35 . 



(7.22) 



Some further numerical results are given in Table V, Again, these parametric 
results should be contrasted with the single' p = 1, v = 2 force-level 
trajectory shown in Figure 3 of [5], 

8 . Final Remarks 

In the previous section above, we have seen how the LCS functions 
allow one to conveniently obtain much valuable information about the model 
(2.1) with power attrition-rate coefficients (3.8) without having to explicitly 
compute the entire force-level trajectories. Previously we were limited to 
computing only force-level trajectories (see [4-5]). With the availability 
of these tabulations of LCS functions (see the Appendix of this report and 
[8]), we can now tell who is going to be annihilated and when this event will 
happen without having to compute the trajectories. Not only did we answer 
questions about the qualitative behavior of the model (e.g, force annihilation) 
for specific values of, for example, initial force levels but also for a 
range of values of the initial force ratio (i.e. parametric analysis of 
model behavior) . 
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Table V. Annihilation of the X Force as a Function 
of the Initial Force Ratio for = 1250 meters 



(x o /y o ) 


X / . N 

t (minutes) 
a 


X , 

r (meters) 
a 


0.333 


10.63 


JL 

1 


0.250 


7.56 


235.9 


0.200 


6.17 


422.8 



t 



t 



max 



9.32 minutes and 



x (r=0) 



1.35. 
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The results of this report may be used for other parametric 
analyses, e.g. parametric dependence of battle outcome on attrition-rate 
coefficients. Thus, the contents of this report allow one to develop 
important insights into the dynamics of combat between two homogeneous 
forces with temporal variations in fire effectiveness. With the avail- 
ability of tabulations of the LCS functions, one can now analyze such 
combat modelled by the power attrition-rate coefficients (3.8) with some- 
what the same facility as he can for the constant-coefficient case and thus 
aid in parametric analyses. For further discussions of the significance 
of such results for military operations research, the redder is directed 
to [6] and [7]. 



34 



REFERENCES 



1. M. ABRAMOWITZ and I. STEGUN (Editors), Handbook of Mathematical 
Functions , National Bureau of Standards Applied Mathematics Series, 

No. 55, Washington, D.C., 1964. 

2. K. KNOPP, Infinite Sequences and Series , Dover, New York, 1956. 

3. F. W. LANCHESTER, "Aircraft in Warfare: The Dawn of the Fourth Arm- 

No. V., The Principle of Concentration," Engineering 98 , 422-423 (1914) 
(reprinted on pp. 2138-2148 of The World of Mathematics , J. Newman 
(Editor), Simon and Schuster, New York, 1956) . 

4. J. G. TAYLOR, "Solving Lanchester-Type Equations for 'Modern Warfare' 
with Variable Cof f icients , " Opns. Res . 22 , 756-770 (1974). 

5. J. G. TAYLOR and G. G. BROWN, "Canonical Methods in the Solution of 
Variable-Coefficient Lanchester-Type Equations of Modern Warfare," 

Opns. Res . 24, 44-69 (1976). 

6. J. G. TAYLOR and G. G. BROWN, "Further Canonical Methods in the 
Solution of Variable-Coefficient Lanchester-Type Equations of Modern 
Warfare," NPS55-77-27, Naval Postgraduate School, Monterey, California, 
June 1977. 

7. J. G. TAYLOR and C. COMSTOCK, "Force-Annihilation Conditions for 
Variable-Coefficient Lanchester-Type Equations of Modern Warfare," 
Naval Res. Log. Quart . 24, 349-371 (1977). 

8. J. G. TAYLOR and G. G. BROWN , "A Table of Lanchester-Clif f ord- 
Schlafli Functions," NPS55-77-39, Naval Postgraduate School, Monterey, 
California, October 1977. 



35 



APPENDIX: Tabulation of the LCS Functions F^ (x) , ^(x), ar >d 

a = 1/2, 1/3, 2/3, 1/4, 3/4, 1/5, 2/5, 3/5, 4/5, 3/7, 



T (x) for 
a 

and 4/7. 
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TABLE 11A. Lanchester-Clif f ord-Schlaf li Functions 
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